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Abstract. We investigate the corrections to scaling of the Renyi entropies of a region 
of size i at the end of a semi-infinite one-dimensional system described by a conformal 
field theory when the corrections come from irrelevant boundary operators. The 
corrections from irrelevant bulk operators with scaling dimension x have been studied 
by Cardy and Calabrese (2010), and they found not only the expected corrections of the 
form £ 4 ~ 2x but also unusual corrections that could not have been anticipated by finite- 
size scaling arguments alone. However, for the case of perturbations from irrelevant 
boundary operators we find that the only corrections that can occur to leading order 
are of the form £ 2 ~ 2x >> for boundary operators with scaling dimension Xb < 3/2, and 
i^ 1 when Xb > 3/2. When Xb = 3/2 they are of the form £ _1 log£. A marginally 
irrelevant boundary perturbation will give leading corrections going as (log^) -3 . No 
unusual corrections occur when perturbing with a boundary operator. 
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1. Introduction 

The block entanglement of a quantum system has been found to be a powerful tool for 
characterizing the scaling behavior near a quantum critical point [T]. For a system in a 
pure state and with the Hilbert space partitioned into a direct product H = Ha ®%b 
(with A and B the corresponding two parts of the system), the block entanglement is 
encoded by the von Neumann entropy Sa = — Trp^ log pA of the reduced density matrix 
Pa, with Sa = Sb- The most interesting case is in one dimension. For an infinite system 
with an interval A of length I the asymptotic behavior of the von Neumann entropy is 
given by [2] 

S A ~^Og-+c[ (1) 

near the critical point. Here c is the central charge of the underlying conformal field 
theory. The constant e is an arbitrary cutoff scale, with c[ also being a non-universal 
number. As a way to characterize the full entanglement spectrum one may introduce an 
additional parameter n, with n a positive real number, and define the Renyi entropies 

^^logTrpl, (2) 

with lim^! = Sa- As expected from finite-size scaling theory, the critical scaling 
S^ ~ (c/6)(l + n~ l ) log(£/e) of the Renyi entropies exhibit 0(£ 4-2:r ) corrections [3]. 
Here x > 2 is the scaling dimension of the leading irrelevant operator (with "irrelevant" 
being understood in the sense of the renormalization group). As shown in Ref. j3], there 
can also be unusual n-dependent corrections of 0(£~ 2x ^ n ) and 0{l 2 ~ x ~ x l n \ where, in the 
first case, x may in fact be less than 2, corresponding to a scaling correction produced by 
a relevant operator. These unusual corrections often come with an oscillating prefactor, 
which however vanishes when n — > 1 in all known cases jH [5]. For a semi- infinite 
system, with a conformally invariant boundary condition (CIBC), operators in the bulk 
may produce additional unusual scaling corrections l~ x l n to the Renyi entropies, on top 
of the ordinary 0(£ 2 ~ x ) corrections with x > 2 predicted by finite-size scaling [3]. In 
contrast to the case of an infinite system, the oscillating prefactor that multiplies the 
leading unusual ^"^-correction does not vanish in the limit n — > 1. This unexpected 
feature was first observed in numerical work in Ref. [6] , and recently derived analytically 
for the case of the XX-chain with open boundary conditions [7]. 

In this article we inquire about the scaling corrections to the critical Renyi entropies 
of a semi-infinite one-dimensional system which are generated by irrelevant boundary 
operators. Recall that boundary operators arise in the operator product expansion 
(OPE) of a chiral operator with its mirror image across the boundary. More precisely, 
given a boundary conformal field theory (BCFT) defined on the complex half-plane 
{z = t + iy | y > 0} with a CIBC at y — 0, the OPE of a chiral operator (j)(r, y) with 
its mirror image 0(r, — y) reads [HI E] 

0(r,y)0(r,-y)~^^p^-^(r), y^ 0. (3) 
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Here x^ is the scaling dimension of 0, and <j)j are boundary operators of dimension Xj. 
Nonzero values of the expansion coefficients C^j select those boundary operators which 
are consistent with the particular CIBC imposed at y = 0. Knowing the boundary 
operator content associated with a system allows for a complete characterization of 
its boundary critical behavior, i.e. those terms in the critical scaling of observables 
contributed by the presence of the boundary. For a quantum theory, where r is a 
Euclidean time, this allows for identifying the long-time (a.k.a. low-energy) asymptotic 
critical behavior of the system close to the boundary. BCFT has a manifold of 
applications, spanning from open-string theory (D branes) [11] to the study of quantum 
quenches [12]. A particularly important class of applications is that of a quantum 
impurity interacting with an electron liquid, where at low energies the impurity can 
be traded for a CIBC at the site of the impurity [10]. The increase of the block 
entanglement at quantum criticality due to the presence of the impurity is a universal 
number (boundary entropy) which characterizes the type of boundary critical behavior. 
However, for a finite block there will always be additive corrections to the boundary 
entropy coming from irrelevant bulk and boundary operators. These corrections are 
expected to reveal features about quantum impurity phenomena which are otherwise 
difficult to access, the extent and character of the enigmatic "screening cloud" being a 
case in point [T3j . 



2. Scaling corrections from bulk operators: a brief review 

Consider a one-dimensional system with a boundary at y = that is described by a 
BCFT. Let subsystem A be the region < y < i and B the rest of the system, y > i. 
As shown by Calabrese and Cardy [21 [I], Tr p n A (which enters the definition of the Renyi 
entropies in Eq. ([2])) can be viewed as a path integral Zn n on an n-sheeted Riemann 
surface lZ n with a boundary, and with proper normalization. Then 

— j— — l°g -j^- = - y~ — - nF ) > (4) 
1 — n Z n 1 — n 

where F = — log Z is the free energy, Fj^ = — log Zn n and (3 is the inverse 
temperature. For an unperturbed BCFT, Zn n /Z n can be calculated as a one-point 
function on the half-plane C + of a twist field $ n inserted at the branch point z = i£ [2] 

7 /9f\ - c ( n - 1 / n )/ 12 

Jr = <<M^)>c+ = c n , (5) 

where e is the short-distance cutoff. This leads to 

4 n) = ^(l + ^)logy + ^ + <, (6) 

where s A = log g A is the boundary entropy and c' n are non-universal constants [21 |T%] . 
The well-known result for the critical block entanglement on a semi-infinite line, 

S A = C -\og- + s A + c>, (7) 
o e 
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is obtained simply by letting n — > 1 + in Eq. (((3]). 

The corrections to this scaling behavior from irrelevant bulk operators was recently 
studied by Cardy and Calabrese [3]. Such a perturbation with a bulk operator $(z) 
having a scaling dimension x > 2 gives an action 

S = Scft + X [ d 2 z<$>(z), (8) 



where A is a coupling constant. The corrections to the free energies are given by the 
perturbation series 



^ = ^-£^r / d V" / d 2 ^^)--'^))^, (9) 



N=l 

over the Riemann surface lZ n , and 

F = F CFT_J2t^l f dW"/ d 2 ^($( Wl )---$(^)) C+ ,(10) 

over the ordinary complex half-plane C + . When the boundary conditions are such that 
($(w))c+ 7^ they found that the first-order correction to Fn n —nF is of the form £ 2 ~ x . 
However, for n > x/(x — 2) they also found the appearance of the unusual n-dependent 
correction £~ x ^ n , that comes from the singularity at the branch point. For an infinite 
system without a boundary, the corrections take the forms £ 4 ~ 2x and £~ 2x l n . 

3. Scaling corrections from irrelevant boundary operators 

We now wish to study the case when the perturbations come from irrelevant operators 
on the boundary. In doing so, we can follow the same procedure as Cardy and Calabrese 
in Ref . [3] . However, when perturbing with a boundary operator the surface integral of 
the perturbing field in the action (jSJ) will be replaced by a line integral on the boundary. 
As we shall see, this will prevent the appearance of unusual n-dependent corrections. 
In fact, this is anticipated since n-dependent exponents only arise from the region at 
the branch point, which is located away from the boundary. Nevertheless, there still 
are results for the boundary case that do not follow from standard finite-size scaling 
analysis. 

Thus, consider a BCFT on the upper half y > of the complex plane z = r + iy, 
so that r is the boundary coordinate at y = 0. The n-sheeted Riemann surface 1Z n 
is then obtained by sewing together n copies of this half-plane along < y < £ at 
r = 0. To evaluate the correlation functions on lZ n for a chiral operator <p(z) with 
scaling dimension Xb, we need to use the transformation property 



w*o =n 



N dz 



((f)(wi) ■ ■ -<f){w N ))c+ , (11) 



dw 

W — Wj 

where the map z w from TZ n to the upper half-plane C + is given by 



( z-U \ V" _ 1 
\z+\t) 1 
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This gives 



dz 
dw 



-Anl 



(w-i\ n 
\w+i) 



|2 - 



(13) 



:i+^ 2 ) [(^) n -i] L 

Naturally, the mapping f lT2|) takes the boundary of lZ n to the boundary of C + . Since w 
is real on the boundary we see from ( {TBI that |dz/du>| _:ri > is analytic on the boundary, as 
the only singularity is at w = i, i.e. when z is at the branch point z = \£. In particular, 
note that the point \z\ — > oo gives a divergence in |dz/dw| which only means that 
\dz/dw\~ Xb ->• 0. 

Now we can use this to study the scaling corrections of Sjf oc (F-ji n — nF) when 
adding a boundary perturbation, 

S = S CFT + X J dr0 b (r), (14) 

where </>& is an irrelevant operator with scaling dimension > 1 on the boundary y = 0. 

We will assume the boundary conditions to be such that (4>b{T)) = 0. This is natural 
if we demand conformal boundary conditions. Then the first-order correction vanishes. 
An important exception is when the perturbing boundary operator is the stress-energy 
tensor, a case that was treated in Ref. pJi|. Since this operator has a non- vanishing 
expectation value on lZ n it will give rise to a first-order correction to \ which was 
found to have the form £~ l . We therefore consider the second-order corrections to F-ji n 
and F, denoted S 2 F-ji n and S 2 F respectively. They are given by 

A* 
'2/3 

and 

X 2 

= --o 

2(3 

respectively, where r[, t' 2 are boundary coordinates on lZ n and T\ = w(r[), r 2 = w(t£) 
are boundary coordinates on C + . From Eq. ( II 1ft . we get 



6 2 F 



dr[ / dr' 2 (MOMO) 



Kn 



dn / dr 2 (0fe(ri)0 6 (r 2 ))c+ 



(15) 



(16) 



S 2 F 



n n - 
2/3 



fdr> 


dz 


-x b 


dz 


-x b 


dw 


z=r[ 


dw 


z=t' 2 





dz 


l—Xb 


dz 


l—Xb 


1 dri J dr 2 


dw 






W=Tl 


dw 


W=T2 



(0 6 (ri)0 6 (r 2 )) c+ . (17) 



We can now use the fact that d^F^ —n5 2 F only depends on the ratio £/e, where e is the 
short-distance cutoff of the theory, to extract its -^-dependence. Since the action (1T4"|) is 
dimensionless, the coupling constant A goes as A ~ e Xb ~ 1 . Thus A 2 ~ e 2x b~ 2 ^ anc i since 
dz/dw oc £ the integral in Eq. (I17p includes an overall factor of (£/e) 2 ~ 2xb . However, 
there can also appear powers of £/e coming from the need to regularize divergences in 
the integrals. 
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In order to compare the two integrals in S 2 F-ji n — nd 2 F it is convenient to rewrite 
n5 2 F on the same form as 5 2 Fn . 



n5 2 F 



:7l 



" 2/3 
2/3 

h n 1 dri 



dn 
n I dri 



dr 9 



1 



d|n - r 2 | 



2.x,., 



Ti - r 2 



2x b 



A-2x b 

- 2x h 



Rewriting this as an integral over the boundary of 1Z n , one gets 



n5 2 F 



dr[ 



A 2 

2/3 



e l-2x 6 

1 - 2x fe 

dz 



dw; 



,1-221!, 



1 - 2x b 



18) 



and then going back to writing this as a double integral over T\ and r 2 gives 



nc, 2 F 



A 2 

2/3 



dn 



dr 9 



ri-T2|>e/|(dz/diu) m= 



d^ 



dw 



2-2x h 



W = T± 



1 



Ti - r 2 



2xh 



As e — > 0, we have 

\(dz/dw) w=Tl \ 



T\ — r 2 > e <^ 



To I > e 



(19) 
(20) 



so that 6~ 2 Fti„ — n5 2 F can be written as a single integral 



6 2 F 



nS 2 F = - 4t / dn / dr- 



2/3 



I dz 
1 d«> 


l—x b 

W=T\ 


dz 
dui 


1— x h 

W = T2 




dz 
dw 


2-2x b 

W=Tl 




Ti 


- r 2 


2x b 



A 2 

4/3 



dri / dr 2 



dz I l ~ x b 



dz \ l ~ x b 



dw 



W = T\ 



dw 



W=T2 



\2.ih 



with the cutoff \(dz/dw] 



W = Tl 



(21) 

l-x b 



IS 



|Ti - r 2 

|ti — r 2 | > e. It follows from Eq. (fT3|) that jdz/dio 
analytic everywhere except at w — ±i. In Ref. j3], where the integrals are over C + , this 
singularity at the branch point was an important ingredient in the analysis. But when 
we now consider a perturbing operator on the boundary, the only divergence in the 
integrand in (12 ip comes when Ti = r 2 . Expanding |dz/d-u;| 1 ~ Xb = f(w) around w = r 2 , 
gives 

(/ / (r 2 )(Ti-r 2 ) + ir(r 2 )(ri-r 2 ) 2 + ...) 2 



5 2 F nn - n5 2 F = — 



A 2 

4/3 
A 2 
4/3 



dn 
dn 



dr 2 



|Ti - r 2 

\2-2x b 



2x b 



dr 2 [(f(r 2 )) 2 |Ti-T 2 | 

+ f (r 2 )f (r 2 )(n - r 2 )|n - r 2 | 2 " 2 ^ + ...]. (22) 

From this it follows that the leading divergence of the double integral in ( 12 lj) goes as 
e 3-2 ^ 6 , i.e. it converges when x b < 3/2. Then no regularization is needed, and the only 
^-dependence comes from dz/dw oc I. Thus, when x b < 3/2 



8 2 F nn - n8 2 F ~ (£/e 



,2-2x b 



(23) 
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and consequently the leading corrections S 2 S^ to the Renyi entropies are of the form 
PS$> ~ f- 2x " . (24) 

On the other hand, when Xb > 3/2, the cutoff in the integral ( 12~TT) must be kept, so that 

5 2 F Un - n5 2 F ~ (£/e) 2 - 2Xb (i/e) 2 ^~ 3 = (l/e)" 1 , (25) 

and then S 2 S^ ~ t~ x for all xj, > 3/2. Note that this is of the same form as the 
first-order correction from the stress-energy tensor. 

When Xb = 3/2, it follows from Eq. (|22"1) that the integral (l2Tj) diverges 
logarithmically, hence 

6 2 S% ] ^rHogi. (26) 

These deviations from the ~ £ 2 ~ 2x b behavior of the correction when Xb > 3/2 goes 
beyond what would be expected from standard finite-size scaling arguments. 

4. Scaling corrections from marginal boundary operators 

When perturbing with a marginal boundary operator one cannot simply put x\, — 1 in 
Eq. (12"T1) and conclude that the second-order corrections to scaling of the Renyi entropies 
vanish, since |dz/dw;| diverges when \z\ — > oo. However, it can be checked that there 
is no need to regularize the integral because of this. We therefore conclude that the 
second-order corrections will be I- independent when Xb — 1. 

Instead of going to the higher-order integrals in the perturbation series of Fn n — nF 
to find the leading ^-dependence of the corrections we will make use of the ^-theorem [16J, 
analogously to how Cardy and Calabrese [3] use the c-theorem in the marginal bulk case. 

The boundary entropy s A of Eq. (J6j) is governed by the "gradient formula" of Friedan 
and Konechny [17] which in our case takes the simple form 

ds A 

ik = ~ m ' (27) 

where is the renormalization-group beta function given by [TS] 

- /3(A) = 4^ = (1 - - TrbX 2 + 0(X 3 ) . (28) 
at 

In the marginally irrelevant case, i.e. with Xb = 1 and X/b > 0, this gives 

^jj = -7T&A 2 + G(A 3 ) , (29) 
with the asymptotic large-£ solution given by 

m ~ Wjry (30) 

Now, as Eq. ( 12?]) becomes ds A /dX = —nbX 2 + 0(X 3 ) when Xb = 1, we have 



s A = const. - ^ A 3 + G(A 4 ) ~ log g A - ~ 3(p/ . • (31) 

3 3n 2 b 2 log {t/ej 



b . q „/./!> . A 1 



Thus the leading correction to the Renyi entropies goes as (log£) 



-3 
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5. Conclusions 

In) 

The result that the leading second-order corrections to the Renyi entropies S A are 
of the form £ 2 ~ 2x b when perturbing with an irrelevant boundary operator with scaling 
dimension Xb < 3/2 holds regardless of the value of n. This is the result anticipated 
from finite-size scaling. When xj, > 3/2 the leading corrections will be of the same 
form as those from the stress-energy tensor, i.e. i~ l . When x& = 3/2 there is also 
a multiplicative logarithmic contribution to the leading correction which then goes as 
£ _1 log£. A marginally irrelevant boundary perturbation gives a correction ~ (log£) -3 . 
Thus there are no unusual n-dependent corrections to scaling of the Renyi entropies 
from boundary operators, as opposed to bulk perturbations where unusual corrections 
to scaling can occur. These unusual corrections originate from the part of the surface 
integral where the bulk operator approaches the branch point created by the Riemann 
surface construction [3]. However, when the perturbing field is on the boundary it never 
comes close to this singularity. 

In Ref. [15] it was found that the first-order correction from the stress-energy tensor 
on the boundary is ~ i^ 1 . It should be noted that this operator is generically present, 
and gives a correction to the entanglement entropy of the same form as that from a 
boundary operator with scaling dimension Xf, > 3/2. 

We also note that the leading corrections of the form £ 2 ~ 2x b from a boundary 
perturbation are similar to what one can get when perturbing with a bulk operator 
in the presence of a boundary. In Ref. [3] it was found that this can give corrections 
of the form £ 2 ~ x , where x is the bulk scaling dimension, but also unusual corrections 
which can dominate. This is therefore a very different situation compared to having the 
perturbing field on the boundary. 
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